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Abstract 

An atomistic approach allowing an accurate and efficient treatment of depolarizing energy and 
field in any low-dimensional ferroelectric structure is developed. Application of this approach 
demonstrates the limits of the widely used continuum model (even) for simple test cases. Moreover, 
implementation of this approach within a first-principles-based model reveals an unusual phase 
transition - from a state exhibiting a spontaneous polarization to a phase associated with a toroid 
moment of polarization - in a ferroelectric nanodot for a critical value of the depolarizing field. 
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Ferroelectric nanostructures (FEN) are of increasing technological and fundamental in- 
terest because of the need in miniaturization of devices, as well as, the appearance of new 
phenomena (see, e.g., Ref. □QDDDQ and references therein). Unscreened polarization- 
induced charges at the surfaces of FEN generate a depolarizing field that is responsible for 
striking properties. Examples are the existence of a critical thickness below which no ferro- 
electricity can appear and the observation and prediction of laminar stripe nanodomains 
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2J, as well as the formation of polarization vortex p, |6| . Interestingly, and despite its huge 
importance, we are not aware of any model being able to exactly calculate the depolarizing 
field and energy in any low-dimensional ferroelectric. For instance, the widely used contin- 
uum model (1) neglects the atomistic nature of materials, (2) is technically applicable only 
in the limit of large enough systems and (3) can not predict the depolarizing energy/field in 
the realistic cases of inhomogeneously polarized samples. 

In this Letter we i) demonstrate that it is possible to derive a scheme allowing the exact 
atomistic computation of the depolarizing energy and field in any low-dimensional FEN; ii) 
use this scheme to check the accuracy of the continuum model for some simple test cases; 
iii) report an unusual phase transition between two different kinds of order parameter in a 
ferroelectric nanodot that is driven by the depolarizing field. 

To calculate the depolarizing energy in low-dimensional ferroelectrics, one first needs 
to realize that a system under perfect open-circuit (OC) electrical boundary conditions 
exhibits a maximum depolarizing field (if the polarization lies along a non-periodic direction), 
while ideal short-circuit (SC) electrical boundary conditions leads to a complete screening 
of charges at the FEN surfaces that fully annihilates any depolarizing field. As a result, the 
depolarizing energy and field experienced by the FEN should involve a difference between the 
dipole-dipole interactions associated with these two extreme electrical boundary conditions. 
We shall write energy of the dipole-dipole interaction in any system in the form 

a/3,ij 

where D = 3, 2, 1 stands for a system periodic in 3, 2 and 1 directions, respectively; D = 
corresponds to non-periodic systems and the sum runs over the atomic sites i and j that 
differ from each other and belong to a supercell (to be denoted by S) mimicking the system. 
Such a supercell is infinitely repeated along the periodic directions, if any. For instance, 
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thin films are modeled by supercells that are repeated in two dimensions while the direction 
associated with the growth direction of the film is non-periodic. For dots, the supercell 
is not repeated. V is the volume of the supercell, p(r*j) the dipole moment at the site i, 
a = x,y,z denotes the Cartesian components. The quantity Q^ S ' D ^ depends on both the 
chosen supercell (S) and the periodicity of the system (D). 

It is straightforward to prove that the elements of the Q matrix for systems periodic in 
three [jj, two (x- and y), one (z) directions 20] and non-periodic systems are given by: 
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where G are the reciprocal lattice vectors associated with the S supercell, A is the Ewald 
ri 

parameter |8[ (which is assumed to be large enough), 5ij is the Kronicker symbol, A is the 
supercell area, p^ and are the projections of (vector connecting atomic sites i andj) on 
the {x,y} plane and z-axis, respectively. T is the incomplete Gamma function and erfc is the 
complementary error function, a is the supercell length in ^-direction; K n are the modified 
Bessel's functions. (Note, that contributions from z,-- n ' )l 



p, L j = in Q ( a p A should be excluded, and that the prime in the right side of Q^ij indicates 
that the term n = has to be excluded when i — j. 

The dipolar interactions described by equations (JT]) and (J2J) correspond to ideal OC condi- 
tions since no charge screening is taken into account in their derivation. The next question 
to be addressed is what is the dipole-dipole energy in FEN under perfect SC conditions. 
Such energy is simply the one described by D = 3 Q-matrix since infinitely extended (i.e, 
bulk) systems do not experience any macroscopic depolarizing field 9|. One can also be 
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convinced in the correctness of the above statement by applying the so-called image method 
that produces polarization pattern identical to the one in the bulk from the polarization 
pattern in the nanostructure. 

The £jZ} (maximum) depolarizing energy per volume in any FEN can now easily be 
calculated as the difference in dipole-dipole energies between perfect OC and SC conditions, 
that is: 

£<Lp = ^dip ~ ^dip = 2Y [Qa/3,ij ~ Q < af3,ij]Pa(. r i)Pf3( r j) (3) 

a/3,ij 

where the sum over i and j run over the sites of the chosen supercell of the FEN, and 
where Q^pij are given by Eq.(J2J). Equation Q is, to the best of our knowledge, the first 
proposed form allowing an atomistic and exact computation of depolarizing energy in any 
FEN with any dipole distribution. (Note that such form can also be applied to calculate 
demagnetization energy in low-dimensional magnetic systems). 

We first apply our approach to compute £j££ in some test cases. Here, we limit ourselves 
to systems adopting a simple cubic structure in which each atomic site has a local and 
equal-in- magnitude dipole moment. Note, that all results to be reported here do not depend 
on the size used for the periodic direction(s) of the S supercell. 

Homogeneous dipole distribution: let's first investigate FEN exhibiting the same local 
dipole p at any atomic site. We shall present our results in the form of = 'jS^ p cont \ 
where is obtained from Eq.(J3j), while S^ con ^ is the depolarizing energy predicted by 
the continuum model. 7 is thus a "correcting" coefficient that provides a measure of the 
continuum approach accuracy. 

For (001) ultra-thin films homogeneously polarized along the out-of-plane (z) direction, 
the continuum model predicts that = 2ttP 2 , independently of the film thickness, 

where P is the polarization. On the other hand, the use of Eq Q results in 7 =1.017, 1.010, 
1.007 and 1.006 for ultra-thin films of 3, 5, 7 and 9 layers, respectively. In other words, 
our atomistic approach reveals that the depolarizing energy is slightly larger than the one 
predicted by the continuum model and increases as the number of film layers decreases. 
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Such findings are consistent with those of Refs.|10llllUl2l|. To understand them, we rewrite 
Eq (JHJ) in the case of a homogeneous dipole pattern as follows: 
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Fig^, shows the "depolarizing" factors < -Fif (I) > averaged over all the j sites belonging 
to a given (001) layer (that is indexed by I) for our films, as a function of a layer position 
inside the film. Comparison with its continuum predicted (and /-independent) value of 4ir 
is also given. FigG^ clearly reveals that the deviation of the continuum model from our 
atomistic results is confined to the surface layers - as also found in Ref. 12] - and that this 

(2) 

deviation is an underestimation. This explains why we numerically found that the £ d J p 
increases with respect to the continuum prediction as the film becomes thinner, since the 
ratio of surface layers over total layers increases as the film thickness decreases. Fig^, also 
shows that for the surface layers < >=< Fyy > are negative - as discussed in Ref. [ll] 

(2) (2) 

- and that 1 + (< F£/ > - < F\V >)/4vr = -0.0393, which is in perfect agreement with 
the calculation of the so-called surface anisotropy in Ref. . 

We next consider different wires of square cross sections, that are periodic along the 
z-direction and homogeneously polarized along the x-axis. According to the continuum 
approach, such wires should have a depolarizing energy £j£f° nt ^ = nP 2 , independently of 
the wire thickness. We numerically found via Eq. (JHJ) that 7 =1.017, 1.010, 1.007 and 1.006 
for wires of 3, 5, 7 and 9 shells, respectively (see inset of FigGJ) for definition of shells). 
Like in the films, the continuum model underestimates the depolarizing energy and this 
underestimation becomes larger in magnitude as the nanostructure shrinks in size. Fig^o 
shows that the continuum model fails to reproduce the averaged depolarizing factor for the 
surface shell but exactly agrees with our atomistic results for all the inner shells. 

We now turn our attention to a cubic dot homogeneously polarized along the z-direction. 
Unlike the previous two cases, our atomistic approach gives a depolarizing energy that is 
not only independent on the dot size but also exactly agrees with the continuum approach 
(that is, £fj p = £f e p° nt ^ = 2ttP 2 /3). Such a surprising result is caused by the sum rule [ill 
for depolarization factors ( that is, F xx + F yy + F zz is the constant given by the continuum 
model) that we numerically found to be valid in every layer (shell) of any system investigated 
so far (i.e. films, wires and dots). Moreover, cubic dots exhibit x, y and z directions that 
are symmetrically equivalent (which is not the case in wires and thin films). As a result, 
< Fxx\l) >=< Fyy (I) >=< FzF\l) > in dots. Because of the sum rule, each of these 
factors is equal to the continuum prediction of 4tc/3 for any shell and for any size, and the 
continuum model predicts the right depolarizing energy. 

Inhomogeneous dipole distribution: let's now investigate 2D, ID and 0D FEN exhibiting 
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stripe domains with the dipoles assumed to be homogeneous inside each domain (with a p 
magnitude) and perpendicular to a periodic direction, if any (see the inset of Fig|2J). The 
period of stripe domains is denoted as d. The £~P} energy calculated from Eq J3J) for a 
film having a thickness L = 10 atomic layers is shown as a function of d/L in Fig|21 along 
with depolarization energies calculated in the continuum approach for the two limiting cases 
d « L (i.e, Sf e 'p° nt) = 1.7P 2 d/L 13, Q) and d » L Q. One can see that these two 
limiting cases can reproduce rather well the energy derived from Eq Q for d/L < 1 and 
d/L > 2, respectively. Furthermore we numerically found that our depolarization energy 
can be parametrized as following 

4ep = Ml - Cl e- d ^ L ) + (nn - c )(l - e~ d ^ L )}P 2 (5) 

where c = 2.568, c x = 1.024, c 2 = 14.118, c 3 = 1.831, n = 2. 

FigJ21 also reports our results for two cases for which we are not aware of any continuum 
predictions, namely stripe domains in infinite wires and cubic dots. Note that the stripe 
direction is along the wire y periodic direction and that the finite size of 0D systems implies 
that d/L has a maximum value of 0.5 in cubic dots. We present here results for a wire of 10 
shells (10x10 atomic sites for cross section) and for a cubic dot of 10 shells (10x10x10 atomic 
sites) with d ranging from 1 to 5 atomic layers. Two features seen in Fig|2]are particularly 
striking. First of all, the stripe domains have less depolarizing energy in a wire than in a thin 
film for the same d/L, with this difference becoming more pronounced as d/L increases. The 
parameters of equation (jSJ) for the wire are c = 2.208, c\ = 1.000, c 2 = 1.196, c 3 = 7.398 
and n = 1. Secondly, for the case of a cubic dot, the dependence of depolarizing energy on 
d/L is linear and given by S^J p = 2.5P 2 d/L. 

We finally take advantage of our formalism to reveal, from realistic first-principles-based 
calculations, properties of a free standing Pb(Zr 0A Ti £)O3 (PZT) cubic dot of 48 A lateral 
size for different electrical boundary conditions. The total energy of the system used in 
Monte Carlo simulations is: 

£ Heff(p(ri), Vi, r],ai) + (3^2 < E dep > p(r^) (6) 

i 

where Sjjeff is the (first-principles-derived Effective Hamiltonian) energy for PZT which 
is dependent on the p(r^) local dipoles at site i of the dot, the Vj inhomogeneous strain 
related variables, the rj homogeneous Btrain tense, a^d on ft e atomi c co^ationQ . 



The dipole-dipole interactions in this H e ff are given by the Q-matrix with D = of Eq|21 
The second term of equation (JHJ) mimics a screening of the (maximum) depolarizing field, 
with the magnitude of this screening being controlled by the (3 coefficient. (3=1 and 
(3 = corresponds to ideal SC and OC electrical boundary conditions, respectively, while 
a value of (3 in-between corresponds to more realistic electrical situation 3j. < Erf ep >= 
— (l/ATcoo) ^2j(d£fo~°^ / dp(rj)) is the depolarizing field inside the dot, while N and eoo are 
the total number of sites of the dot and the dielectric constant of PZT, respectively. 
is practically calculated from Eq. (|3*|). 

FigOk and b show the resulting macroscopic dipole moment and the macroscopic toroid 
moment of polarization (i.e., the supercell average of the cross product between position and 
dipole moment |6j]), respectively, as a function of (3. One can clearly see that for situations 
close to SC, the dot exhibits a macroscopic polarization, with a cross section of the local 
dipole pattern being given in the inset of FigO^. On the other hand a dot with electrical 
boundary conditions close to OC has a non- vanishing toroid moment[Q\, with a cross-section 
of the corresponding dipole pattern being displayed in the inset of FigOb- Moreover, FigOl 
clearly reveals that, at a critical value of the depolarization field, the system undergoes 
an unusual phase transition between a state characterized by one kind of order parameter 
(toroid moment of polarization) to a state associated with another kind of order parameter 
(polarization). In other words, no coexistence between these two order parameters occurs. 

In summary, we have derived an atomistic, simple, general and efficient approach to 
calculate the depolarizing energy and field in any low-dimensional ferroelectric structure. 
The application of this method reveals - and explains - the limits of the continuum model, 
and also results in the discovery of an unusual phase transition in a ferroelectric dot for 
some critical value of the residual depolarizing field. 
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1-0600 and by DOE grant DE-FG02-05ER46188. 
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FIG. 1: Depolarizing factors < Fxx\l) >, < Fyy \l) > and < F£\l) > obtained with our 

atomistic approach and normalized to Att (that is the prediction of F zz in the continuum model) in 

the case of (001) films (part (a)) and 2ir (prediction of F xx in the continuum model) in the case of 

wires periodic along the z-axis (part (b)) as a function of the layer or shell index /. The different 

shells of a wire are shown in the inset of part (b). The most inner layer or shell is indexed by 0. 
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FIG. 2: Normalized depolarization energy as a function of d/L for quasi 0D, ID and 2D systems 
with stripe domains along with continuum model prediction for d « L ^| and d » L 3] for 
D=2 systems. The inset shows the schematic representation of the chosen polarization pattern in 
(001) films, or wires periodic along y or cubic dots. 
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FIG. 3: a) Dipole moments and b) toroid moment of polarization G in a PZT nanodot as a 
function of the screening coefficient (3. Insets of the parts a) and b) show the polarization pattern 
for (5 = 1 (SC conditions) and (3 = (OC conditions), respectively. 
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